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ABSTRACT 


This  paper  presents  a  method  for  calculating  the  probability  that 
a  certain  minimum  number  of  a  set  of  I CBM's  can  be  destroyed  by  a  given 
number  of  interceptor  satellites. 
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PROBABILITY  CONSIDERATIONS  ON  DESTROYING  ICBM'S 
WITH  INTERCEPTOR  SATELLITES 

The  problem,  as  defined  by  Dr.  Fritz  Hoehndorf ,  may  be  formu¬ 
lated  as  follows: 

Let  m  and  n  be  the  numbers  of  interceptors  and  ICBM's,  respec¬ 
tively.  Each  of  the  m  interceptors  is  supposed  to  attack  exactly 
one  of  the  n  ICBM's.  One  ICBM  may  be  attacked  by  one  or  more  in¬ 
terceptors.  The  probability  that  one  interceptor  destroys  the 
attacked  ICBM  is  p.  The  probability  that  the  I  interceptor 
attacks  the  k  ICBM  shall  be  independent  of  i,  k. 

What  is  the  probability  that  at  least  n  -  ^  ICBM’s  are  de¬ 
stroyed  or,  in  other  words,  at  most  ^  ICBM's  escape  destruction? 

Let  us  first  answer  the  question:  What  is  the  probability  of 
destroying  exactly  J  ICBM's?  Let  mx  be  the  number  of  interceptors 
attacking  the  first  ICBM,  let  m2  be  the  number  of  interceptors  at¬ 
tacking  the  second  ICBM,  and  finally  let  mQ  be  the  number  of  inter- 
ceptors  attacking  the  n  ICBM.  It  is  obvious  that  the  number  of 
different  possibilities  of  attributing  mx>  m2,  ...  interceptors 
to  the  first,  second,  ...  nth  ICBM  is 

_ al  .  - 
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Since  nm  is  the  total  number  of  possible  distributions  of  m  "ele¬ 
ments”  among  n  "boxes,"  the  distribution  n^,  nig,  ...  1%  has  the 
probability 


ml 


nm  n^I  n^!  ...  m^l 

The  probabilities  that  the  first,  second,  ...  ICEM's  escape 

nig  mj 

destruction  are  q  ,  q  ,  ...  q  ,  respectively,  where  ?»1  -  p. 

The  probabilities  that  the  (j+l)st,  ...  n^  ICEM's  are  destroyed 
nu.,  bl. 

are  1  -  q  ...  1  -  q  -  ,  respectively. 

Hence  the  configuration;  first  3  ICEM's  not  destroyed,  remain¬ 
ing  n  -  j  ICEM's  destroyed,  itil  ICEM  attacked  by  m^  interceptors 
(i  *»  1,  ...  n),  has  the  probability 

\  - 7  ^  ...  (i  -  &)  ...  Cx  -  <xV) 

n“  S  mgj  ...  m^l  \  /  \  / 

Since  there  are  0)  possibilities  of  selecting  a  set  of  j  escapes 
out  of  the  n  ICEM's  we  obtain 

(1) 

where  is  the  probability  that  exactly  J  ICEM's  escape  destruction 
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and  where  the  sum  has  to  he  extended  over  all  configurations 
m^  nig,  ...  nijj  so  that 

mx  +  nig  +  •••  =  m 

The  stun  in  formula  (l)  may  he  written  in  a  more  convenient  form. 
We  substitute  in  equation  (l)  ' 

%  +  nfe  +  •••  mj  -  v 

Since  each  term  of  the  sum  in  equation  (l)  for  which  mj+k  **  0 
(k  2  l)  vanishes,  we  have  to  consider  only  those  configurations 
for  which 

♦•••%■»-  0^  ♦  •••  »J>  -  B  -  V  *  n  -  ) 
or 


v  S  m  -  n  +  j 


The  set  of  the  v  elements  n^,  . ..  mj  can  he  selected  in  ^  ™  ^ 
ways  and,  there  are  possihilities  of  attributing  v  interceptors 
to  3  ICBC's.  The  number  of  possihilities  of  attributing  the  re¬ 
maining  m  -  v  interceptors  to  the  remaining  n  -  J  ICBC's  is 


(m  -  v)l 
mJ+x*  •••  “n1 
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where 


/  V  f  md+i\  f  ®n\ 

any  mj+k  *  0  (k  *  1).  If  we  now  replace  (  1  -  q  )  ...  (  1  -  q  ) 


by  its 


My  under  the  boundary  condition  (3)  we  see  that 


S  *  My 


I 


(m  -  v)  l 


(5) 


where  only  those  configurations  have  to  be  considered  for  which 


"j+n  *  1 

This  leads  to  the  following  problem:  How  many  possibilities  are 
there  to  distribute  P  "elements"  to  a  "boxes"  (P  >  o)  so  that  each 
box  contains  at  least  one  of  the  P  elements? 

In  order  to  answer  this  question  we  shall  derive  a  recursion 
formula  as  follows:  Let  *p(a,  j)  be  the  number  of  possible  dis¬ 
tributions  bo  that  the  first  j  boxes  contain  at  least  one  element. 

It  is  evident  that  Sp (a,  j+l)  equals  ep(cr,  j)  minus  that  number 
of  distributions  for  which  the  first  J  boxes  contain  at  least  one 
element  and  for  which  the  (j+l)8t  box  is  empty.  This  number  is 
♦p(c-l,  j).  Therefore: 

♦p (o,  J+l)  -  ♦?(<*,  J)  -  ♦pfff-l,  J)  (6) 

We  now  have 

♦p (o,  0)  -  aP  (7) 
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From  equations  (6)  and  (7)  we  conclude  "by  induction 

a 

♦p (or,  <0-  0**  -  c(c-l)P+  ^  (oHB)**  -  •••«  (-l)V  (o-v)P  (8) 

v=0 

or 

a 

♦p(»>  «)  -  £  <-Dv  (  "  )  Vp  (9) 

VaG 

where  Op (a,  a)  means  the  number  of  possibilities  to  distribute  P 
"elements"  to  a  "boxes"  so  that  each  box  contains  at  least  one  ele¬ 
ment.  It  is  desirable  to  express  *p(o,  a)  in  a  closed  form.  For 
this  purpose  we  write 

£  <-Dv  (  l  )  «VI-  (.*  -  D°  (10) 

v=*0 


Differentiating  (10)  P  times  with  respect  to  x  we  obtain 


♦p  (o,  a) 


(•*  -  l)ff 


(11) 


x  =  0 


which  yields  Op  (a,  <r)  in  the  desired  closed  form. 
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Equation  (11 )  shows  that  *p(o,  a)  may  be  written  as  a  contour 
integral  in  the  Gaussian  plane,  extended  over  any  closed  contour 
encircling  the  origin  once 


♦p(c,  a) 


(12) 


In  various  ways  we  can  obtain  useful  estimates  for  $p (c,  a),  for 
instance  choosing  the  unit  circle  about  the  origin  as  contour. 
However,  we  shall  not  dwell  upon  this  method  any  longer,  because 
the  fundamental  formula  (2)  can  be  transformed  into  an  interesting 
and  much  simpler  form,  handy  above  all  for  computations  with  a 
digital  computer. 

With  regard  to  practical  calculations  the  greatest  difficulty 
in  evaluating  formula  (2)  is  caused  by  the  factorial  terras  and, 
above  all,  by  decomposing  m  into  a  large  number  of  summands.  How 
this  can  be  avoided  shall  now  be  investigated. 

The  second  sum  on  the  right-hand  side  of  equation  (2)  has  the 
structure 


(13) 


where 


®1  ^  st  °  8 
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Eliminating  the  parentheses  in  (13)  we  have  ^  ^products  of  the 

1 

S,  So  6/y  f  t  \ 

kind  q  q  ...  q  w  .  All  these  f  )  products  contribute  the 


same  value 


_  /  .  \0  V  si  si  8a 

S  »  (_l)  )  — - —  q  ...  4 

^  8«^  •  •  •  •  8^« 


(i*0 


to  the  sum  in  expression  (13)  if  we  first  keep  s1,  s2,  ... 
fixed.  We  substitute  sx  +  sa  +  ...  sa  »  s  and  write  (l4)  in  the 
form 


Keeping  first  Blt  sz,  ...  sa  fixed  and  summing  over  s^,  ...  st  , 
we  obtain 


(-l)a  (t  -  a)8"®  q8 


(16) 
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Summing  over  sx,  ...  aa  ,  and  keeping  s  fixed  during  this  sum¬ 
mation,  ve  obtain 

(-1)“  (t  -  c.)8'8  45  (  !  )  a8"  (17) 

or 

S  -  (-l)a  (t  -  a  p)8  (18) 

Equation  (l&)  now  leads  to  the  Important  result 


s  3  z 


(19) 


CM) 


Equation  (19)  finally  enables  us  to  give  an  expression  for  Pj 
much  more  suitable  for  practical  calculations: 


m-n+J  n-j 

■. -Oil (:)(<■)' z '->“(■?)(- -»r  <*> 


v-0 


CM) 
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Let  us  now  represent  Pj  as  a  contour  integral.  For  this  purpose 
we  apply  to  the  second  sum  in  equation  (20)  a  method  analogous  to 
that  which  led  to  equation  (11)#  From  the  relation 


I  ™°t3) 


(n-j-ap)x  (n- J )x 


(l  - 


we  conclude  by  differentiating  (m  -  v)  times  with  respect  to  x,  and 
subsequently  substituting  x  «  0 


or 


I 


(n-d) 


(m-v)i 


.n-v+1 


dz  (25) 


where  the  contour  Integral  has  to  be  extended  over  a  closed  contour 
about  the  origin  in  the  Gaussian  plane  of  the  complex  variable  z. 


10 


Substituting  equation  (23)  in  equation  (20)  we  obtain  after  some 


elementary  operations 


m-n+j 


f  . jl/  y^ 

J  \«V  nm  2i*  ^  ^  vl 

v=0 


(> 


-pz 


T'J 


m+1 


dz  (24) 


The  sum 


m-n+j 

I 


v-0 


(jqz)v 

vl 


(25) 


can  rigorously  be  replaced  by  e^2  because  powers  of  jqz  higher  than 
m-n+j  would  not  yield  any  residue  different  from  zero  in  the  integral 
of  equation  (24).  Hence  it  follows  that  Pj  can  be  represented  in  the 
form 


It  is  appealing  to  check  equation  (26)  by  proving  that  indeed 

n 

£  Pj  -  1  (27) 

J-0 
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Since 


j-0 

ve  obtain  from  equation  (26) 


(28) 


Since 


/ 


„nz  » 

e  4z  m  2ix  5L 


equation  (27)  has  been  verified. 

Equation  (26)  can  be  written  in  the  form 


Pj  -  (  “  )  .  A.  .  SL  f 

°  d  >  n*  2ix  « 


(eZ  -  (>  .  XJ 


.B+l 


dz 


(29) 
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It  suggests  eliminating  the  parentheses  in  equation  (26)  in  order 
to  calculate  the  residues.  After  some  elementary  transformations  we 
obtain: 


For  practical  computations  it  is  always  advisable  to  comprehend  two 
succeeding  terms  In  the  sum  of  equation  (30): 


--U*. 

n-3P 


n-J-u 

U+l 


— 2—  T} 

n-Jp-uP  /  J 


(3D 


^  h  _  a-J-u  e“  a-jp-ptp 

n“JP  /  V  n+1 


This  approximation  is  sufficient  for  all  practical  purposes.  Finally 
equation  (30)  assumes  the  form; 


i 


o*0 


n-J 

2a 


(32) 
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for  an  odd  (n-j).  For  an  even  (n-j)  the  upper  limit  of  the  sum  equals 

(n  -  j) 

2 

Professor  Harry  Carver  suggested  that  ve  calculate  the  factorial 
moments  of  our  distribution  function,  remarking  that  the  general  ex¬ 
pressions  would  probably  become  very  simple. 

The  Vth  factorial  moment  is  defined  as  follows: 

H(v)  »  £  J(J  -  1)  {i  -  V  +  1)  Pj  (33) 

J=0 

or  using  the  notation: 

JO-1)  ...  (j  -  v  +  l)  -  j(v) 

we  have: 

“(v)  *  I  j(v)  PJ  <34) 

>0 

In  order  to  calculate  these  moments  we  make  use  of  formula  (26)  which  we 
shall  write  in  a  slightly  modified  form: 


lb 


From  equation  (35)  we  conclude': 


(>  .  3 

\n 

il  v  j<v'  D 

2*1 

L3  \iJ 

iepz-lT'!  iz 


J- 0 


(36) 


The  sum  in  the  integrand  can  easily  be  evaluated  as  follows:  We 
substitute, 

epz  -  1  -  x'1 

Since; 


e  (;)(•*■•>)"•  1  (;)*•■  (‘*-T 

J=0  j=0 


(37) 


we  obtain  by  deriving  equation  (37)  V  times  with  respect  to  x  and 
multiplying  the  result  by  xV: 


Ij<v)  cy*”(v)*v 

J“0 


(38) 


or 


f  J(V)  (  J  )-  (fZ  -  l)"J  "  n(v)  ePZ(n”V)  (ePZ  -  l)’a  (39) 

J-o 
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1  ,1 
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We  now  substitute  equation  (39)  in  equation  (36): 


(v) 


^Cv) 


ml  .  nvv/  r 
2tn  nm  J 


(n  -  yp)z 


m+i 


dz 


(40) 


Calculating  the  residue  of  the  integrand  of  equation  (40)  we 
finally  obtain: 


«(v)-*(v)0-?7 


m 


Our  next  aim  is  to  give  an  approximation  of  our  distribution 
function  good  to  any  desired  degree  of  accuracy,  taking  advantage 
of  the  simple  expressions  for  the  factorial  moments. 

We  are  particularly  interested  in  the  case  where  n  is  very 
large  i.e.,  n  >  200.  In  this  case  it  is  justified  to  replace  the 
summation  sign  in  equation  (34)  by  an  integral  sign: 

n 

“<*)  "  /  J<V)  p(j) 

,3=0 

We  now  perform  the  similarity  transformation: 

J-f  (*  +  1)  (*3> 
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/which  transform*  the,  interval  (0,  n)  into  the  interval  (-1,  l),  and 
obtain: 


(kb) 

•V 

where  p(x)  dx  is  the  probability  that  x  lies  in  the  interval 
(x,  x  +  dx). 


Now  we  take  advantage  of  the  simple  expressions  for  the  factorial 
moments  in  order  to  expand  P(x)  into  a  series  of  Legendre  polynomials. 
For  this  purpose  we  eliminate  the  parentheses  in  equation  (44): 


/i  (l  +  x)  P(x)  dx  -  |  s3L(v"1)  J'  ^  (l*  xf1  p(x) 


where  si^v_1)  means  the  i"**  elementary  symmetric  function  of  the 
first  v  -  1  integers. 
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For  practical  applications  ve  shall  list  the  first  four  symmetric 
functions  sx^n^,  s2^;  83^,  sa^  : 

8l(*o * 


s  (n)  -  n(n  +  il  (n  zJ^LSM  +  _2) 

2  "  24 

(46) 

S3(n)  =  alia  ±iJLila 

48 


S4(n)  =  n(n  +  1)  (n  -  l)  (n  -  2)  (n  -  3)  (15n3  +  15na  -  lQn  -  8) 

27  •  3*  *  5 


Equation  (45)  suggests  expanding  the  expression'  (l  +  x)n  into  a  sum  of 
Legendre  polynomials.  As  will  he  shorn  in  the  Appendix,  we  find: 


a  ♦*)”■=  £  ita  pt 

i«0 


(*7) 


tth 


where  is  the  i  Legendre  polynomial  normalized  so  that: 


fn 

-i 


dX  a  1 


(48) 
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(21  +  1)1  |^(n+l)n  -  (i+l)ij  ^n(n-l)  -  (l+l)ij  ...^(i+2)  (i+l)  -  (i+l)  i 
Let  us  now  expand  p(x)  into  a  series  of  Legendre  polynomials: 


P(x)  «  Y  (x)  (5°) 

6 

We  nov  substitute  equations  (4l),  (47),  and  (50)  in  equation  (45): 


iv  i 


a.  - 


(v-i) 


n 


i=0 


L 

1-0 


iv-i 


(51) 


+  £ 


V— 2 


I 


iv-2  i 


1=0 


Equation  (5l)  enables  us  to  calculate  easily  the  coefficients  a^j. 
successively,  and  is  above  all  suitable  for  large  numbers  n  and  m. 


Although  equation  (51 )  apparently  is  complicated,  it  permits 
computation  of  the  coefficients  a^  without  solving  a  system  of  linear 
equations.  The  numerical  evaluation  of  equation  (51)  can  be  performed 
by  means  of  a  desk  calculator. 


^frsswsswsts 


APPENDIX 

In  order  to  prove  formula  (49) ,  page  19,  we  start  from  the  dif¬ 
ferential  equation  for  Legendre's  polynomials: 


dx 


(x2  -  1)  P^  »  i(i  +  1)  P± 


(1) 


We  now  multiply  equation  (l)  "by  (l  +  x)n  and  integrate  from 
-1  to  +1  : 

+i  i 

i(i  +  1)  J  (1  +  x)n  P±  dx  *  f  §£  (x2  “  1)  P[  (1  +  x)n  dx 


=  -  n  fX  (x2  -  1)  (1  +  x)n~x  p;  dx  -  -  n  /%-  1)  (1  +  x)n  P.' 


dx 


(2) 


i  i 

n  J  (1  +  x)n  Pt  dx  +  n2  J  (x  -  l)  (l  +  x)n~x  T±  dx 

-l  -l 

(n  +  n2)  J  (i  +  x)n  dx  -  2a2  J  (l  +  x)11”1  dx 


From  equation  (2)  we  derive  the  recursion  formula: 


Iin  = 


2n2 


n2  +  n  -  i2  +  i 


"in-i 


(3) 


21 
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I 

I 


! 


k 


! 
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We  apply  formula  (3)  successively  until  ve  reach  1^,  which  is: 


hi  "  f  (!  +  x)1  dx  -  f  x1  Pi 


dx 


(4) 


The  integral  (4)  is  veil  known  and  may  he  taken  from  any  textbook*-, 
and  from  it  ve  then  obtain,  after  some  elementary  transformations, 
equation  (49),  page  19. 
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*W.  Magnus  and  F.  Oberhettinger,  Formeln  und  Satze  fur  die  speziellen 
Funktionen  der  mathematischen  Physik. 
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